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We study electronic structures of quasi-two-dimensional finite electron systems in high magnetic 
fields. The solutions in the fractional quantum Hall regime are interpreted as quantum liquids of 
electrons and vortices. The ground states are classified according to the number of vortices inside 
the electron droplet. The theory predicts observable effects due to vortex formation in the chemical 
potentials and magnetization of electron droplets. We compare the transitions in the theory to those 
found in electron transport experiments on a quantum dot device and find significant correspondence. 
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Studies of finite two-dimensional electron systems in 
high magnetic fields often get their inspiration from the 
remarkable physics of the quantum Hall effect in two- 
dimensional electron gas y). Considerable experimen- 
tal and theoretical work has been carried out in par- 
ticular on quantum dots, two-dimensional droplets of 
electrons in the interface region of semiconductor het- 
crostructures 0, 0, 0. 0|- Electron transport measure- 
ments of quantum dots have revealed a rich variety of 
transitions associated with charge redistribution within 
the electron droplet in high magnetic fields 0. This 
has lead to development and study of theoretical mod- 
els which could account for the microscopic origin of 
these phenomena 0, EJ El- Recently a proposal 
has been put forth that the observed phenomena may 
be caused by emergence of a special state of electronic 
matter, a quantum liquid of strongly correlated electrons 
and vortices EI 0|- In the present work we 

use this approach to classify the ground states of the 
two-dimensional electron droplets in quantum dots and 
compare transitions to those found in electron transport 
measurements. Our results give also general insight into 
the internal structure of the many-body wave function of 
finite two-dimensional electron systems. 

Vortices, rotational flow of currents or matter with a 
characteristic cavity at the center, can be found in vari- 
ous natural phenomena where particles have been set to 
rotate around a common axis. Descriptions of vortices 
exist since antiquity [lfij . In two-dimensional quantum 
systems a vortex can be defined in analogue with classi- 
cal vortices as a zero in the wave function associated with 
a phase change of integer multiple of 2ir for each path 
enclosing this zero. In quantum dots the rotation of elec- 
trons is induced by external magnetic field and vortices 
may form if this rotation is strong. Vortices are caused 
by quantization of the magnetic flux through the electron 
droplet and they give rise to rotating currents of charge 
around density zeros inside the electron droplet. Vor- 
tices create charge deficiency inside the electron droplet 
which manifests itself as an increase of the dot area (jjla]- 
This theory has been based on electronic structure stud- 



ies of quantum dots [ill El El as well as on theoretical 
analogies with bosonic systems of rotating Bose-Einstein 
condensates 0. Vortices in electron droplets are not 
necessarily bound to electrons as approximated by the 
Laughlin wave functions. The analysis of the internal 
structure of the many-body wave function suggests intro- 
duction of a more general framework of interacting sys- 
tem of electrons and off-electron vortex quasi-particles, 
where vortex formation is driven by interactions 0, E| • 
We use an effective-mass approximation in the xy- 
plane to model the physics of quasi-two-dimensional sys- 
tem of trapped electrons in vertical quantum dot devices. 
The electron-electron interaction is approximated with 
Coulomb potential. The Hamiltonian is then 
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where N is the number of electrons, V c is the external 
confining potential, m* is the effective mass of electrons 
moving in semiconductor medium, e is the dielectric con- 
stant, and A is the vector potential of the homogeneous 
magnetic field which is oriented perpendicular to the xy- 
plane. In the subsequent discussion the external poten- 
tial is chosen to be parabolic V c (r) = \m*J^r 2 . 

In this Letter we use both the mean-field spin- 
density-functional theory (SDFT) and variational quan- 
tum Monte Carlo (VMC) to calculate the electronic 
structure in the quantum Hall regime. We use the SDFT 
in conjunction with local spin density approximation 
(LSDA) with a smooth correlation functiona l Il7l. For 
details of the implementation we refer to Refs. |lllll8l |. In 
high magnetic fields the electron droplet is spin-polarized 
and a stable stable structure called the maximum den- 
sity droplet (MDD) forms El- If i s a finite size precur- 
sor of the integer v = 1 quantum Hall state. Evidence 
for the MDD formation has been reported in the experi- 
ments 0, . When the magnetic field is increased in the 
measurements the MDD breaks down into a lower density 
droplet in the fractional quantum Hall (FQH) regime. 

Theory predicts that in parabolically confined quan- 
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FIG. 1: Mean-field SDFT phase diagrams of the ground state 
of parabolically confined electron droplets in high magnetic 
fields. The ground states are classified according to the num- 
ber of vortices inside the electron droplet (gray-scale). The 
confining potential is hu>o = 5 meV in a) and 7.67-/V -1 / 7 meV 
in b). In the upper right hand corner of the right diagram the 
number of vortices is greater than 8. 

turn dots the ground states in the beyond-MDD (FQH) 
regime occur only at certain "magic" angular momen- 
tum values |2*oj . As a result the angular momentum as a 
function of external magnetic field shows a characteristic 
staircase structure. In the MDD breakdown the angu- 
lar momentum increase with respect to the MDD state 
AL = L — L(MDD) strongly depends on the number of 
electrons in the system |l8|. For N < 12 the electron in 
the center is moved to the outer edge giving AL = N 
and a vortex hole emerges at the center of the dot. For 
N > 12 a vortex emerges at a finite distance from the 
center This change in the breakdown mechanism 
means that vortices in large electron systems tend not 
to localize. However, if the symmetry of the external 
potential is broken localized vortices may form in the 
particle and current densities |l4| . High angular momen- 
tum states correspond to multi-vortex configurations in 
the FQH regime. Since they are beyond reach for exact 
diagonalization techniques for N > 10 we use the SDFT 
to analyse the electronic structure of these states. 

In the experimental realizations of quantum dots the 
area of the dot has been found to increase with the gate 
voltage suggestin g th at the electron density in the dot 
remains constant |2l| . In zero magnetic field this implies 
a confining potential that scales as Hujq ~ iV -1 / 4 |22j. 
However, the magnetic field exerts additional squeezing 
effect on the electrons which is counteracted by interac- 
tions. We have found that approximately constant elec- 
tron density in the calculations is obtained by a potential 
scaling Huq ~ N~ 1/>7 . Fig.njshows phase diagrams of the 
ground states in the SDFT. The ground states are clas- 
sified according to number of vortices inside the electron 
droplet by using a conditional single-determinant wave 
functions constructed from Kohn-Sham orbitals 0. 

For qualitative understanding of the structure of the 
phase diagram we interpret the solutions in the FQH 
domain as a quantum liquid of electrons and hole-like 
vortex quasi-particles. The MDD state assigns one Pauli 
vortex at each electron position. As the MDD recon- 



structs the number of vortices in the system increases by 
one and subsequent transitions involve emergence of more 
off-electron vortices. Since vortices carry magnetic flux 
quanta <I>o the MDD state breaks down approximatively 
when the magnetic flux $ = BA through the MDD of 
area A exceeds (A + l)$o and an additional off-electron 
vortex emerges in the electron droplet. The compact- 
ness of the MDD state and the observed constancy of 
the electron density n with respect to A gives A = nA 
and B = (1 + 1/A)$ n for the MDD reconstruction. The 
upper boundary of the MDD is therefore approximatively 
constant for high TV in accord with the experiments |6j. 
In the addition of one vortex in the MDD the relative 
increase of the dot area is approximately 1/N. This in- 
crease of the dot area via vortex formation counteracts 
the squeezing effect of the magnetic field. Assuming now 
constant n for the beyond-MDD states the change in the 
B required for addition of subsequent off-electron vortices 
in the droplet is approximately AB = ^^n/N . The cal- 
culations show that widths of the different vortex phases 
are to a good accuracy inversely proportional to A which 
gives credence to this simplified picture (see Fig.^ 3 ). The 
above reasoning implies also that constancy of the elec- 
tron density in the droplet leads to narrowing of the MDD 
window for high electron numbers. The lower boundary 
of the MDD is determined either from polarization of 
the electron droplet or flipping of one vortex from paral- 
lel to anti-parallel orientation with respect to the mag- 
netic field. The latter condition involves a constant flux 
change of 2$o with respect to the MDD state. This leads 
to narrowing of the MDD window because the area of 
the droplet increases with A in a totally polarized sys- 
tem. However, for realistic Zeeman-coupling strengths 
the electron droplet loses total spin-polarization in the 
low-field limit before antivortices emerge. This shrinks 
the MDD window further. 

We now compare these theoretical results with ex- 
perimental data from electron transport measurements. 
Oosterkamp and coworkers measured electron transport 
through a vertical quantum dot device in the quantum 
Hall regime Up to this date these experiments give 
the best available electron transport data for a single-dot 
device. Fig. shows chemical potentials in these exper- 
iments for electron numbers N — 12 to 39, the MDD 
window boundaries in the SDFT as well as the transi- 
tions associated with emergence of off-electron vortices. 
The gate voltage dependence of the external confinement 
in the quantum dot device is taken into account by using 
confining potential hu Q = 5.70A" 1 / 7 meV in the SDFT. 
For material parameters we use m* — 0.067, e = 12.4, 
and effective gyromagnetic ratio g* — —0.44. Despite 
the simple form of the inter-electron potential in our the- 
oretical model the transitions in the beyond-MDD do- 
main fit well into the experimental data. This can be un- 
derstood from the fact that the transitions are induced 
by the magnetic flux quantization through the relatively 
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compact electron droplet. In Fig. [3] the chemical poten- 
tial of the 24-electron quantum dot is calculated with 
the SDFT, the VMC, and the exact diagonalization in 
the lowest Landau level (LLL). The results are compared 
to the experimental data in Ref. |(J . The correspondence 
between the theory and experiments is good and differ- 
ent Quantum Hall regimes can be identified by comparing 
the two sets of data. We find also that the theoretical 
results are consistent with the observed narrowing of the 
MDD window from about IT at N = 20 to 0.5 T at 
N = 39 (see Fig.EJ- In the SDFT the partially polarized 
states in the spin-flip region have a MDD-like compact 
structure for both orientations of the z-component of the 
spin. Therefore the states before the MDD around 6 T 
for high N in Fig. could be partially polarized states 
and the MDD window may be smaller than that identi- 
fied in Ref. [6j . In the theory the increase in the dot area 
at the MDD breakdown is approximately l/N=3.3% at 
N = 30. This can be contrasted to around 10% reported 
in the experiments @|. However, uncertainty in the ex- 
perimental result is high and there exists no analysis of 
a possible N dependence [26|. 

The electron transport data shows features which 
may be due to correlation effects beyond our mean-field 
SDFT. These may include the transition associated with 
the open triangle in Fig. and fluctuations in the data 
between the second and third triangle. In addition, the 
LLL theory of Ref. suggest that for small particle num- 
bers (< 100), the first ground state after MDD would 
have partial spin polarization. However, higher Landau 
levels might have effect on this [2^. In actual quan- 
tum dot realizations, effects due to finite thickness of the 
electron gas, screening of the interaction potential, image 
charges, and non-parabolic terms in the external poten- 
tial may also cause deviations from our results |24|. 

In the electron transport experiments the slope of 
the chemical potential is the magnetization difference 
M(N,B)-M(N-1,B), where M= -BE/3B. Since the 
magnetic field couples to the angular momentum L, there 
is a jump in the magnetization at transitions associated 
with change in L. In the MDD state the angular momen- 
tum increases with N as N — > N+l. Therefore the slopes 
of the chemical potentials in the MDD domain show typ- 
ical decreasing pattern as N increases at fixed B. SDFT 
calculations indicate that magnetization per electron in 
the FQH regime is close to the effective Bohr magneton 
H* B = m e /m*n B = 0.864 meV/T (see Fig. EJ). This is 
a consequence of emergence of additional vortex quasi- 
particles in the electron droplet as the magnetic field in- 
creases. This behaviour gives rise to positive slopes in the 
chemical potentials in the beyond-MDD domain which is 
in accord with experimental data (see e.g. slopes of the 
chemical potentials at B = 9 T in Fig. [2J ■ Detail plots of 
the magnetization for N = 5 and N = 24 in the mean- 
field theory (Figs. 0Jd and reveal also oscillations in 
the droplet magnetization as the number of vortices in- 




FIG. 2: Current peaks in the electron transport experiments 
and transitions in the SDFT (red lines). The dashed lines 
denote the MDD boundaries. The right dashed line and the 
dotted lines correspond to transitions associated with emer- 
gence of off-electron vortices one-by-one inside the electron 
droplet. The experimental data is from Fig. 2b in Ref. |(|. 
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FIG. 3: Chemical potential of the 24-electron quantum dot 
calculated with the SDFT and compared to the experiments 
from Ref. and VMC and LLL exact diagonalization results 
in the vicinity of the MDD window. Noise in the experimen- 
tal data has been reduced by using a Gaussian filter. The 
numbers between the arrows indicate the total spin. The 
roman numbers between the triangles indicate the number 
of vortices inside the electron droplet in the fractional quan- 
tum Hall regime. Charge-density-wave solutions in the SDFT 
(dotted line) have been discarded as unphysical. 

creases one- by-one 14]. The sawtooth- like oscillations 
in lower magnetic fields in Fig. ^Ji are manifestations of 
the de Haas-van Alphen effect. The overall behaviour is 
similar to that found in direct measurements of magneti- 
zation of dot mesas |27| . This method could also provide 
a way to detect oscillations in the FQH regime. 

To conclude, we have performed theoretical calcula- 
tions for finite electron droplets in the Quantum Hall 
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FIG. 4: a) Ground state magnetization per electron calculated 
with the SDFT. The confining potential is 7.677V" 1/7 meV. 
The beyond-MDD states are characterized by a plateau re- 
gion in this plot, where magnetization per electron is close to 
\x" B . The dashed and the dotted lines correspond to the same 
transitions as in Fig. |2] b) and c) Detail plots of the oscilla- 
tions in the magnetization of the 5-electron and 24-electron 
droplets, respectively. The finite-size precursor of the v = 1/3 
quantum Hall state is identified by using conditional wave 
functions [Tl|. The confinement strength is 5 meV in a) and 
3.62 meV in b). The oscillations in the plot a) are smoothed 
out due to use of finite temperature. 



regime. Our model theory predicts emergence of a quan- 
tum liquid of electrons and off-electron vortices in high 
magnetic fields. The pattern of transitions found within 
the theory is consistent with experimental data from elec- 
tron transport measurements. However, for greater qual- 
itative understanding of the phenomena we call for more 
accurate electron transport experiments in the quantum 
Hall regime, more accurate computational methods, and 
more realistic modeling of quantum dot systems in the 
case of many-electron systems. It has been suggested 
that addition of cusps to smooth LSDA functionals incor- 
porates more FQH correlations 28] . Another approach 
towards testing of the validity of theoretical predictions 
could be a direct visualization of the electron density in a 
quantum dot. Recent developments with scanned probe 
imaging techniques |29|, |30j could yield a way to image 
localized vortices in electron droplets. 
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